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In this work we present new exact similarity solutions with moving boundaries of the 
Fokker-Planck equation having both time-dependent drift and diffusion coefficients. 
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§1. Introduction 



One of the basic tools which is widely used for studying the effect of fluctuations 
in macroscopic systems is the Fokker-Planck equation (FPE).O This equation has 
found applications in such diverse areas as physics, chemistry, hydrology, biology, 
finance and others. Because of its broad applicability, it is therefore of great interest 
to obtain solutions of the FPE for various physical situations. 

Generally, it is not easy to find analytic solutions of the FPE, except in a few 
simple cases, such as linear drift and constant diffusion coefficients. In most cases, 
one can only solve the equation approximately, or numerically. Most of these meth- 
ods, however, are concerned only with FPEs with time-independent diffusion and 
drift coefficients (for a review of these methods, see eg. Ref. [JJ)). 

Solving the FPEs with time-dependent drift and/or diffusion coefficient is in 
general an even more difficult task. It is therefore not surprising that the number of 
papers on such kind of FPE is far less than that on the FPE with time-independent 
^ coefficients. Some recent works on the FPE with time-dependent diffusion coefficients 

appear in and works involving time-dependent drift coefficients can be 

found in[5]),[6D,|7D. Refs. ED J9]) JTOD consider FPEs with both time-dependent diffusion 
and drift coefficients. The symmetry properties of the one-dimensional FPE with 
arbitrary coefficients of drift and diffusion are investigated in ITU]) . Such properties 
may in some cases allow one to transform the FPE into one with constant coefficients. 

In Ref. 7) we have considered, within the framework of a perturbative approach, 
the similarity solutions of a class of FPEs which have constant diffusion coefficients 
and small time-dependent drift coefficients. The solvability of the FPE with both 
time-dependent drift and diffusion coefficients by means of the similarity method 
is explored in Ref. [TT]) . One advantage of the similarity method is that it allows 
one to reduce the FPE to an ordinary differential equation which may be easier to 
solve, provided that the FPE possess proper scaling property under certain scaling 
transformation of the basic variables. Some interesting exactly solvable cases of such 
FPE on the real line x € (—00,00) and the half lines x £ [0,oo) and x € (— 00, 0] 
were considered in Ref. Hip. These domains admit similarity solutions because their 
boundary points are the fixed points of the scaling transformation considered. This 
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indicates that similarity solutions are not possible for other finite domains. 

It is therefore natural to suspect that similarity solutions of FPE on a finite 
domain may be possible, if its boundary points scale accordingly. We are thus led to 
consider FPE with moving boundaries. The purpose of this note is to present three 
new classes of such FPEs. To the best of our knowledge, FPEs with time-dependent 
coefficients and moving boundaries have not been discussed before. 



§2. Scaling of Fokker-Planck equation 



We first recapitulate the scaling form of the FPE. The general form of the FPE 
in (1 + l)-dimension is 



dW(x,t) 
Ft 
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where W(x, t) is the probability distribution function, D^\x, t) is the drift coefficient 
and D^ 2 \x,t) the diffusion coefficient. The drift coefficient represents the external 
force acting on the particle, while the diffusion coefficient accounts for the effect of 
fluctuation. W(x, t) as a probability distribution function should be normalized, i.e., 
/domain W ( x > t) dx = 1 ioi t > 0. 

We shall be interested in seeking similarity solutions of the FPE, which are 
possible if the FPE possesses certain scaling symmetry. 

Consider the scale transformation 

x = e a x , t = e b t, (2-2) 

where e, a and b are real parameters. Suppose under this transformation, the prob- 
ability density function and the two coefficients scale as 

W(x,t) =e c W(x,t), D^(x,t) = e d D^(x,t), D^(x,t) = e e D^(x,t). (2-3) 

Here c, d and e are also some real parameters. It can be checked that the transformed 
equation in terms of the new variables has the same functional form as eq. (|2-1|) if the 
scaling indices satisfy b = a — d = 2a — e. In this case, the second order FPE can be 
transformed into an ordinary differential equation which may be easier to solve. Such 
reduction is effected through a new independent variable (called similarity variable), 
which is certain combinations of the old independent variables z such that it is 
scaling invariant, i.e., no appearance of parameter scaling transformation is 

performed. Here the similarity variable z is defined by 

x a 
z = — , where a = — and a , b ^ . (2-4) 

For a , b ^ 0, one has a ^ , oo. 

The general scaling form of the probability density function W(x, t) is W(x, t) = 
(t Sl /x Xl )y(z), where Ai and 5\ are two real parameters and y(z) is a scale-invariant 
function under the same transformation (|2-2p . From the assumed scaling behavior 
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of W(x,t) in (|2-3p . we have —c/a = X% — (5%/a). Without loss of generality and for 
clarity of presentation, we take the parameters (Ai,#i) = (0, etc/ a). This gives 

W(x,t) = t<y(z) , (2-5) 

where y(z) is a function of z. The normalization of the distribution function is 

e^y(z) 



domain 



W(x,t)dx 



domain L 



dz = 1 



(2-6) 



For the above relation to hold at all t > 0, the power of t should vanish, and so one 
must have c = —a, and thus 



W{x,t) 



(2-7) 



Similar consideration leads to the following scaling forms of the drift and diffusion 
coefficients 



D^(x,t) = t a - 1 Pl (z) , D^(x,t) = t 2a - 1 p 2 (z) , 

where pi(z) and p 2 (z) are scale invariant functions of z. 
With eqs. (15% I^B and tiffy . the FPE is reduced to 



(2- 



P2 (z) y"{z) + [2p' 2 (z) - Pl (z) + az\ y'(z) + \p'i{z) - p[(z) + aj y(z) = , (2-9) 

where the prime denotes the derivative with respect to z. It is really interesting to 
realize that eq. (|2-9p is exactly integrable. Integrating it once, we get 



P2(z)y'(z) + [p' 2 {z) - pi(z) + az] y(z) = C, 
where C is an integration constant. Solution of eq. (|2-10p is 



(2-10) 



y( z ) 



c' + c 



dz- 



- F dzf(z) ' 



cxp 



dzf(z] 



t( \ _ Pi( z ) - P2( z ) ~ az ( n / n 

f{z) = r~\ ' P2K Z )TQ, 

P2[Z) 



(2-11) 



where C is an integration constant. 

To proceed further, let us consider the boundary conditions of the probability 
density W(x,t) and the associated probability current density J(x,t). From the 
continuity equation 



dt 



dx 



d r. 



J(x, t) = (x, t) W(x, t)-— D (2) (x, t) W(x, t) 

ox L 



we have 



Using eqs. (jlM]) . (p 7 ?) and (j^ . we get 

J(x, t) = r 1 [( Pl (z) - p' 2 (z)) y(z) - p 2 (z)y'(z)] 



(2-12) 
(2-13) 

(2-14) 
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From eq. (|2T0p . we can reduce the above equation to 




(2-15) 



where W(x,t) = t~ a y(z) and z = x/t a have been used in obtaining the second 
expression. 

In this paper we consider boundaries which are impenetrable. Hence we must 
have 



These conditions imply that C = 0, and that J(x,t) is proportional to W(x,t) and 
x. 

With C = 0, the probability density function W(x,t) is given by 



where A is the normalization constant. It is interesting to see that the similarity 
solution of the FPE can be given in such an analytic closed form. Exact similarity 
solutions of the FPE can be obtained as long as p\{z) and P2(z) are such that 
the function f(z) in eq. (|2-17p is an integrable function and the resulted W(x,t) 
is normalizable. Equivalently, for any integrable function f(z) such that W(x,t) is 
normalizable, if one can find a function P2{z) (pi(z) is then determined by f(z) and 
P2(z)), then one obtains an exactly solvable FPE with similarity solution given by 
(|2-17p . Some interesting cases of such FPE on the real line x £ (—00,00) and the 
half lines x £ [0, 00) and x S (—00, 0] were discussed in Ref. [TT]) . 

That the half and whole lines can accommodate similarity solutions which allow 
x to scale is obvious: the points x = and x = ±00 are the fixed points of the scale 
transformation ()2-2|) . Other finite intervals do not allow similarity solutions, unless 
their end points scale accordingly. We are thus led to FPE with moving boundaries. 

In what follows we shall present three classes of exactly solvable FPE of such 
kind. 



We consider a finite domain x\(t) < x < X2(t) with impenetrable moving bound- 
aries at Xk(t) (k = 1,2). We want the transformed FPE in z-space can be exactly 
solvable. So the simplest choice is such that in the z-space the boundary points 
of the corresponding domain are static. This implies Z}~ = Xk(t)/t a (k = 1,2) are 
constants. We will assume this choice below. We note here that the fixed do- 
mains admitting similarity solutions considered in Ref. [TT]) correspond to the choice 
z = x(t)/t a = 0, ±00, which are just the fixed points of the scaling transformation. 

Now let us assume the function f(z) to have the following form in the physical 
domain in the z-space: 




(2-16) 




(2-17) 



§3. Class I: two moving boundaries 



/(*) 



Ol, Cl2 > 0, Z\ < Z < Z2- 



(3-1) 



Z — Z\ 



Z2 - Z 
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This leads to a choice of pi(z) and p 2 (z): 

p 2 {z) = (z- Zl )(z2 - z), (3-2) 
pi(z) = (a - ai - a 2 - 2)z + (ai + l)z 2 + (a 2 + l)zi (3-3) 

for z\ < z < z 2 , and pi(z), p 2 {z) = otherwise. The function y(z) in the physical 
domain is 

y{z) = A{z-z 1 )^{z 2 -z) a \ (3-4) 

Here the normalisation constant A is given by 

A = [(z 2 - Zl ) a ^ +1 B( ai + l,a 2 + (3-5) 

where B(x,y) is the Beta function. 
The probability density function is 

w {x ,t)={ t (#-*i) ai (*-#r> 

l 0, otherwise 
There are three subclasses: 

(i) Zl , z 2 >0 ( Zl , z 2 < 0); (3-7) 

(ii) z x = 0, z 2 > (zi < 0, z 2 = 0); (3-8) 

(iii) < 0, z 2 > 0. (3-9) 

The situations given in the brackets correspond to mirror images of the corresponding 
classes with an appropriate change of parameters. In subclass (ii), z\ = is a fixed 
point of the scale transformation, and can be considered as a special case of Case II 
to be discussed below. 

In Figs. 1 to 3 we show figures for subclass (i) and (iii). For Figs. 1 and 3 we set 
a > 0, while for Fig. 2 we have a < 0. It is seen that for a > 0(< 0), the boundaries 
move away from (toward) the origin (except when the boundary is a fixed point). 
The same pattern is also exhibited in the other classes. 

§4. Class II: One moving boundary with x = a fixed point 

We now consider the case on the positive half-line with x = a fixed point. The 
case with the moving front in the negative half-line is simply the mirror images of 
the situation discussed here with an appropriate change of parameters. 

The function y(z) is given by 

f( z ) = — ^— + 0, ai, a 2 > 0, < z < Z2, P ■ real. (4-1) 

z z 2 - z 

This leads to a choice of pi(z) and p 2 (z): 

p 2 (z) = z{z 2 - z), (4-2) 
pi(z) = -(3z 2 + (a-a 1 -a 2 -2 + /3z 2 )z + (en + l)z 2 (4-3) 
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W(x,t) 




Fig. 1. W(x, t) vs x for Case-I(i) with a = 2, zi = 1, zi = 4, ai = 1 and ci2 = 1/2 for t = 0.3 (dotted 
line), 0.4 (dashed line) and 0.5 (solid line). 

W(x,t) 




\ 

i 
i 
i 

1 2 3 



Fig. 2. W(x,t) vs x for Case-I(i) with a = -2,zi = \,zi = 4, ai = 1/3 and a 2 = 1/2 for t = 1.0 
(dotted line), 1.2 (dashed line) and 1.4 (solid line). 



for < z < Z2, and pi, p2 = otherwise. The function y(z) is 

y{z) = Az ai {z 2 - z) a2 e^ z , (4-4) 

where the normalisation constant A is given by 

A = [^ 1+ ° 2+1 B(ai + l,a 2 + l) iF 1 (ai + l;a 1 + a 2 + 2;/3z 2 )]- 1 , (4-5) 

where \F\{n\ u; z) is the Kummer confluent hypergeometric function. 
The probability density function is 

W(x t) = { £ ^ ~ e ^< ° * x * **** (4-6) 

1 0, otherwise 
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Fig. 3. W(x, t) vs x for Case-I(iii) with a = 2, z\ = —2, Z2 — 4 and ai = a2 = 1 for i = 0.6 (dotted 
line), 0.8 (dashed line) and 1.0 (solid line). 



A representative situation of this class is depicted in Fig. 4. 



W(x,t) 
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0.6 
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Fig. 4. W(x, t) vs x for Case-II with a = 2, Zi = 1, a\ = 1, aa = 1/2 and /3 
line), 0.6 (dashed line) and 0.8 (solid line). 



-1 for t = 0.4 (dotted 



§5. Class III: One moving boundary with x = oo a fixed point 

This last class represents the case on the positive half-line with x = oo a fixed 
point. The situation with x = — oo a fixed point is the mirror image of the situation 
discussed here with an appropriate change of parameters. 
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For this class the function y{z) is taken to be 
a\ ci2 



f3, a\, a2 > 0, z\ < z, f3 > 0. 



Z — Z\ Z 2 

The corresponding choice of pi(z) and P2{z) is: 

P2(z) = (Z - Zi)Z, 

pi(z) = -fiz 2 + (a + a\ + a 2 + 2 + fiz\)z - (a 2 + l)zi 
fo z\ < z, and pi, p2 = otherwise. The y(z) is 

y(z) = A(z-zi) ai z a2 e- /?2 . 
Here the normalisation constant A is given by 

A=[^ ~z l 2 r(ai + l)e-g^ gl Wa.-a, a 1+a2+ i (jgza)]" 1 , 

2 ' 2 

where -T(2i) and W^ )U (z) are the Gamma and Whittaker functions, respectively 
The probability density function is 



(5-1) 



(5-2) 
(5-3) 



(5-4) 



(5-5) 



W(x,t) 



(#) e "^- zit Q <x<oo 
0, otherwise 



(5-6) 



Fig. 5 shows the evolution of W(x, t) for some choice of parameters. 



W(x,t) 
1.2 
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Fig. 5. W(x,t) vs x for Case-Ill with a = 2, 2i = 0.5, ai = 1, a 2 = 1/2 and /3 = 1 for £ = 0.5 
(dotted line), 0.8 (dashed line) and 1.0 (solid line). 



In summary, we have found that similarity solutions of the FPE are possible 
with moving boundaries, provided that the boundaries points scale appropriately. 
Three new classes of exactly solvable FPE's with moving boundaries are presented. 
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